ELULC-10, a 10 m European land use and land cover map using Sentinel and landsat data in Google Earth Engine

Land Use/Land Cover (LULC) maps can be effectively produced by cost-effective and frequent satellite observations. Powerful cloud computing platforms are emerging as a growing trend in the high utilization of freely accessible remotely sensed data for LULC mapping over large-scale regions using big geodata. This study proposes a workflow to generate a 10 m LULC map of Europe with nine classes, ELULC-10, using European Sentinel-1/-2 and Landsat-8 images, as well as the LUCAS reference samples. More than 200 K and 300 K of in situ surveys and images, respectively, were employed as inputs in the Google Earth Engine (GEE) cloud computing platform to perform classification by an object-based segmentation algorithm and an Artificial Neural Network (ANN). A novel ANN-based data preparation was also presented to remove noisy reference samples from the LUCAS dataset. Additionally, the map was improved using several rule-based post-processing steps. The overall accuracy and kappa coefficient of 2021 ELULC-10 were 95.38% and 0.94, respectively. A detailed report of the classification accuracies was also provided, demonstrating an accurate classification of different classes, such as Woodland and Cropland. Furthermore, rule-based post processing improved LULC class identifications when compared with current studies. The workflow could also supply seasonal, yearly, and change maps considering the proposed integration of complex machine learning algorithms and large satellite and survey data.Peer ReviewedPostprint (published version

Ocean remote sensing techniques and applications: a review (Part II)

As discussed in the first part of this review paper, Remote Sensing (RS) systems are great tools to study various oceanographic parameters. Part I of this study described different passive and active RS systems and six applications of RS in ocean studies, including Ocean Surface Wind (OSW), Ocean Surface Current (OSC), Ocean Wave Height (OWH), Sea Level (SL), Ocean Tide (OT), and Ship Detection (SD). In Part II, the remaining nine important applications of RS systems for ocean environments, including Iceberg, Sea Ice (SI), Sea Surface temperature (SST), Ocean Surface Salinity (OSS), Ocean Color (OC), Ocean Chlorophyll (OCh), Ocean Oil Spill (OOS), Underwater Ocean, and Fishery are comprehensively reviewed and discussed. For each application, the applicable RS systems, their advantages and disadvantages, various RS and Machine Learning (ML) techniques, and several case studies are discussed.Peer ReviewedPostprint (published version

Mangrove Ecosystem Mapping Using Sentinel-1 and Sentinel-2 Satellite Images and Random Forest Algorithm in Google Earth Engine

Mangroves are among the most productive ecosystems in existence, with many ecological benefits. Therefore, generating accurate thematic maps from mangrove ecosystems is crucial for protecting, conserving, and reforestation planning for these valuable natural resources. In this paper, Sentinel-1 and Sentinel-2 satellite images were used in synergy to produce a detailed mangrove ecosystem map of the Hara protected area, Qeshm, Iran, at 10 m spatial resolution within the Google Earth Engine (GEE) cloud computing platform. In this regard, 86 Sentinel-1 and 41 Sentinel-2 data, acquired in 2019, were employed to generate seasonal optical and synthetic aperture radar (SAR) features. Afterward, seasonal features were inserted into a pixel-based random forest (RF) classifier, resulting in an accurate mangrove ecosystem map with average overall accuracy (OA) and Kappa coefficient (KC) of 93.23% and 0.92, respectively, wherein all classes (except aerial roots) achieved high producer and user accuracies of over 90%. Furthermore, comprehensive quantitative and qualitative assessments were performed to investigate the robustness of the proposed approach, and the accurate and stable results achieved through cross-validation and consistency checks confirmed its robustness and applicability. It was revealed that seasonal features and the integration of multi-source remote sensing data contributed towards obtaining a more reliable mangrove ecosystem map. The proposed approach relies on a straightforward yet effective workflow for mangrove ecosystem mapping, with a high rate of automation that can be easily implemented for frequent and precise mapping in other parts of the world. Overall, the proposed workflow can further improve the conservation and sustainable management of these valuable natural resources

Linear and Non-Linear Vegetation Trend Analysis throughout Iran Using Two Decades of MODIS NDVI Imagery

Vegetation is the main component of the terrestrial Earth, and it plays an imperative role in carbon cycle regulation and surface water/energy exchange/balance. The coupled effects of climate change and anthropogenic forcing have undoubtfully impacted the vegetation cover in linear/non-linear manners. Considering the essential benefits of vegetation to the environment, it is vital to investigate the vegetation dynamics through spatially and temporally consistent workflows. In this regard, remote sensing, especially Normalized Difference Vegetation Index (NDVI), has offered a reliable data source for vegetation monitoring and trend analysis. In this paper, two decades (2000 to 2020) of Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI datasets (MOD13Q1) were used for vegetation trend analysis throughout Iran. First, the per-pixel annual NDVI dataset was prepared using the Google Earth Engine (GEE) by averaging all available NDVI values within the growing season and was then fed into the PolyTrend algorithm for linear/non-linear trend identification. In total, nearly 14 million pixels (44% of Iran) were subjected to trend analysis, and the results indicated a higher rate of greening than browning across the country. Regarding the trend types, linear was the dominant trend type with 14%, followed by concealed (11%), cubic (8%), and quadratic (2%), while 9% of the vegetation area remained stable (no trend). Both positive and negative directions were observed in all trend types, with the slope magnitudes ranging between −0.048 and 0.047 (NDVI units) per year. Later, precipitation and land cover datasets were employed to further investigate the vegetation dynamics. The correlation coefficient between precipitation and vegetation (NDVI) was 0.54 based on all corresponding observations (n = 1785). The comparison between vegetation and precipitation trends revealed matched trend directions in 60% of cases, suggesting the potential impact of precipitation dynamics on vegetation covers. Further incorporation of land cover data showed that grassland areas experienced significant dynamics with the highest proportion compared to other vegetation land cover types. Moreover, forest and cropland had the highest positive and negative trend direction proportions. Finally, independent (from trend analysis) sources were used to examine the vegetation dynamics (greening/browning) from other perspectives, confirming Iran’s greening process and agreeing with the trend analysis results. It is believed that the results could support achieving Sustainable Development Goals (SDGs) by serving as an initial stage study for establishing conservation and restoration practice

Satellite-Observed Spatial and Temporal Sea Surface Temperature Trends of the Baltic Sea between 1982 and 2021

The Baltic Sea is one of the fastest-warming marginal seas globally, and its temperature rise has adversely affected its physical and biochemical characteristics. In this study, forty years (1982–2021) of sea surface temperature (SST) data from the advanced very high resolution radiometer (AVHRR) were used to investigate spatial and temporal SST variability of the Baltic Sea. To this end, annual maximum and minimum SST stacked series, i.e., time series of stacked layers of satellite data, were generated using high-quality observations acquired at night and were fed to an automatic algorithm to detect linear and non-linear trend patterns. The linear trend pattern was the dominant trend type in both stacked series, while more pixels with non-linear trend patterns were detected when using the annual minimum SST. However, both stacked series showed increases in SST across the Baltic Sea. Annual maximum SST increased by an average of 0.062 ± 0.041 °C per year between 1982 and 2021, while annual minimum SST increased by an average of 0.035 ± 0.017 °C per year over the same period. Averaging annual maximum and minimum trends produces a spatial average of 0.048 ± 0.022 °C rise in SST per year over the last four decades

A Dual Attention Convolutional Neural Network for Crop Classification Using Time-Series Sentinel-2 Imagery

Accurate and timely mapping of crop types and having reliable information about the cultivation pattern/area play a key role in various applications, including food security and sustainable agriculture management. Remote sensing (RS) has extensively been employed for crop type classification. However, accurate mapping of crop types and extents is still a challenge, especially using traditional machine learning methods. Therefore, in this study, a novel framework based on a deep convolutional neural network (CNN) and a dual attention module (DAM) and using Sentinel-2 time-series datasets was proposed to classify crops. A new DAM was implemented to extract informative deep features by taking advantage of both spectral and spatial characteristics of Sentinel-2 datasets. The spectral and spatial attention modules (AMs) were respectively applied to investigate the behavior of crops during the growing season and their neighborhood properties (e.g., textural characteristics and spatial relation to surrounding crops). The proposed network contained two streams: (1) convolution blocks for deep feature extraction and (2) several DAMs, which were employed after each convolution block. The first stream included three multi-scale residual convolution blocks, where the spectral attention blocks were mainly applied to extract deep spectral features. The second stream was built using four multi-scale convolution blocks with a spatial AM. In this study, over 200,000 samples from six different crop types (i.e., alfalfa, broad bean, wheat, barley, canola, and garden) and three non-crop classes (i.e., built-up, barren, and water) were collected to train and validate the proposed framework. The results demonstrated that the proposed method achieved high overall accuracy and a Kappa coefficient of 98.54% and 0.981, respectively. It also outperformed other state-of-the-art classification methods, including RF, XGBOOST, R-CNN, 2D-CNN, 3D-CNN, and CBAM, indicating its high potential to discriminate different crop types

Spatial and Temporal Analysis of Surface Urban Heat Island and Thermal Comfort Using Landsat Satellite Images between 1989 and 2019: A Case Study in Tehran

Mapping and monitoring the spatio-temporal variations of the Surface Urban Heat Island (SUHI) and thermal comfort of metropolitan areas are vital to obtaining the necessary information about the environmental conditions and promoting sustainable cities. As the most populated city of Iran, Tehran has experienced considerable population growth and Land Cover/Land Use (LULC) changes in the last decades, which resulted in several adverse environmental issues. In this study, 68 Landsat-5 and Landsat-8 images, collected from the Google Earth Engine (GEE), were employed to map and monitor the spatio-temporal variations of LULC, SUHI, and thermal comfort of Tehran between 1989 and 2019. In this regard, planar fitting and Gaussian Surface Model (GSM) approaches were employed to map SUHIs and derive the relevant statistical values. Likewise, the thermal comfort of the city was investigated by the Urban Thermal Field Variance Index (UTFVI). The results indicated that the SUHI intensities have generally increased throughout the city by an average value of about 2.02 °C in the past three decades. The most common reasons for this unfavorable increase were the loss of vegetation cover (i.e., 34.72%) and massive urban expansions (i.e., 53.33%). Additionally, the intra-annual investigations in 2019 revealed that summer and winter, with respectively 8.28 °C and 4.37 °C, had the highest and lowest SUHI magnitudes. Furthermore, the decadal UTFVI maps revealed notable thermal comfort degradation of Tehran, by which in 2019, approximately 52.35% of the city was identified as the region with the worst environmental condition, of which 59.94% was related to human residents. Additionally, the relationships between various air pollutants and SUHI intensities were appraised, suggesting positive relationships (i.e., ranging between 0.23 and 0.43) that can be used for establishing possible two-way mitigations strategies. This study provided analyses of spatio-temporal monitoring of SUHI and UTFVI throughout Tehran that urban managers and policymakers can consider for adaption and sustainable development

Performance Evaluation of Six Gridded Precipitation Products throughout Iran Using Ground Observations over the Last Two Decades (2000–2020)

Precipitation, as an important component of the Earth’s water cycle, plays a determinant role in various socio-economic practices. Consequently, having access to high-quality and reliable precipitation datasets is highly demanded. Although Gridded Precipitation Products (GPPs) have been widely employed in different applications, the lack of quantitative assessment of GPPs is a critical concern that should be addressed. This is because the inherent errors in GPPs would propagate into any models in which precipitation values are incorporated, introducing uncertainties into the final results. This paper aims to quantify the capability of six well-known GPPs (TMPA, CHIRPS, PERSIANN, GSMaP, IMERG, and ERA5) at multiple time scales (daily, monthly, and yearly) using in situ observations (over 1.7 million) throughout Iran over the past two decades (2000–2020). Both continuous and categorical metrics were implemented for precipitation intensity and occurrence assessment based on the point-to-pixel comparison approach. Although all metrics did not support the superior performance of any specific GPP, taking all investigations into account, the findings suggested the better performance of the Global Satellite Mapping of Precipitation (GSMaP) in estimating daily precipitation (CC = 0.599, RMSE = 3.48 mm/day, and CSI = 0.454). Based on the obtained continuous metrics, all the GPPs had better performances in dry months, while this did not hold for the categorical metrics. The validation at the station level was also carried out to present the spatial characteristics of errors throughout Iran, indicating higher overestimation/underestimation in regions with higher precipitation rates. The validation analysis over the last two decades illustrated that the GPPs had stable performances, and no improvement was seen, except for the GSMaP, in which its bias error was significantly reduced. The comparisons on monthly and yearly time scales suggested the higher accuracy of monthly and yearly averaged precipitation values than accumulated values. Our study provides valuable guidance to the selection and application of GPPs in Iran and also offers beneficial feedback for further improving these products

Flood Hazard Mapping Using Fuzzy Logic, Analytical Hierarchy Process, and Multi-Source Geospatial Datasets

Iran is among the driest countries in the world, where many natural hazards, such as floods, frequently occur. This study introduces a straightforward flood hazard assessment approach using remote sensing datasets and Geographic Information Systems (GIS) environment in an area located in the western part of Iran. Multiple GIS and remote sensing datasets, including Digital Elevation Model (DEM), slope, rainfall, distance from the main rivers, Topographic Wetness Index (TWI), Land Use/Land Cover (LULC) maps, soil type map, Normalized Difference Vegetation Index (NDVI), and erosion rate were initially produced. Then, all datasets were converted into fuzzy values using a linear fuzzy membership function. Subsequently, the Analytical Hierarchy Process (AHP) technique was applied to determine the weight of each dataset, and the relevant weight values were then multiplied to fuzzy values. Finally, all the processed parameters were integrated using a fuzzy analysis to produce the flood hazard map with five classes of susceptible zones. The bi-temporal Sentinel-1 Synthetic Aperture Radar (SAR) images, acquired before and on the day of the flood event, were used to evaluate the accuracy of the produced flood hazard map. The results indicated that 95.16% of the actual flooded areas were classified as very high and high flood hazard classes, demonstrating the high potential of this approach for flood hazard mapping